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Abstract. We present preliminary results of the generalized Prin- 
cipal Component Analysis (PCA) of light curves of 82 magnetic 
chemically peculiar (further mCP) stars applied to 54 thousand in- 
dividual photometric observations in the uvby and H p colours taken 
from the "On-line database of photometric observations of mCP 
! stars" . 



We show that every of the observed light curves (LCs) can be, 
with a sufficient accuracy, represented using five parameters of a 
harmonic polynomial of the second order, and that the third order 
and higher harmonics reflect only a noise. We found that a prevailing 
majority of the uvbyH v LCs can be satisfactorily well represented by 
a linear combination of a constant term and one or two (extraordi- 
narily three) basic, mutually orthonormal functions, which reduces 
the number of parameters necessarily needed to describe the set of 
the LCs of a particular star. 

While the shape of an individual LC depends on the individual 
distribution of photometric spots on the surface of a rotating star, its 
amplitude depends on the mechanism of variability. To describe the 
amplitude of LCs we introduce a robust quantity effective amplitude. 
Applying the PCA to the set of the effective amplitudes of the LCs 
in all five colours we revealed at least three different sources of the 
light variability showing different wavelength dependence. The am- 
plitudes of all the LCs harbor a component monotonously decreas- 
ing with increasing wavelengths. The second component reaches a 
remarkable extremum in the v colour, and is in antiphase relative to 
the remaining four colours. It occurs particularly at cooler mCP-s 
having a high degree of the chemical peculiarity. The third compo- 
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nent reflects a remarkable diversity of the light curves in the u colour 
found in a few mCP stars. 

Key words: light variability - light curves - chemically peculiar 
stars - origin of the light variations 



1 Introduction 

The light variability of the magnetic chemically peculiar (mCP) stars, ranging a few hun- 
dredths, seldom a few tenths of magnitude in amplitude, has been known and investigated 
for decades. These strictly periodic light variations, frequently accompanied by synchronous 
variations of spectral lines and magnetic field, are commonly explained due to existence of 
hypothetical, vast and persistent photometric spots on the surface of rigidly rotating stars. 
The properties of these photometric spots and their relation to the known spectroscopic spots 
and the geometry of the magnetic field have not been explored satisfactorily up to now. 

We aim at a deep analysis of the rotationally modulated photometric variations of mCP 
stars, their astrophysical explanation based both on the complex study of selected well ob- 
served mCP stars and the study of the photometric behaviour of as large as possible sample 
of the stars. 



2 The data 

Prior to start the investigation we established a publicly accessible "On-line database of 
photometric observations of mCP stars" (http://dumbell.physics.muni.cz/pilot, hereafter the 
Database) as today containing more than 120 thousand observations of more then 120 mCP 
stars made in 40 photometric colours. One can find a more detailed information on this public 



database in the paper Mikulasek et al.(2007) 



We selected 82 stars having observations in u, v, b, y and Hipparcos H p colours. The sam- 
ple contains one He-strong (HD 37776), 10 He- weak, 37 Si and 34 "cool" SrCrRE stars. Here 
we study altogether 408 LCs in five photometric colours consisting of 54 131 individual ob- 
servations with errors of ~ 0.003 mag. 

The linear light elements of all the 82 stars were determined by means of our own method 
based on the generalized PCA and a robust regression, based on the complete photometric 
material contained in the Database. The fastest rotator in our sample is HD 164429 with the 
rotational period P = 0.519 d, the slowest one is HD 188041 with P = 224 d. The median is 
P = 3.8 d. 

An evident trend that the cooler mCP stars rotate slower on average, can be seen on Fig. 
[TJ The undereddened B- V indices were taken from the Hipparcos photometry. 



3 Description of mCP light curves 

The rotational LCs of magnetic CP stars are strictly periodic what means they can be de- 
scribed to advantage by harmonic polynomials depending on the rotational phase (p. Further, 
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Figure 1: Dependence of the rotational period on the mean Hipparcos B-V index. 

it follows from the experience (INorth 19841) that the rotational LCs are smooth to such ex- 
tent that the polynomial of the 2nd degree is quite sufficient. The terms of the higher degrees 
mostly reflect noise so they can be neglected. Thus 

m(c,(p) = a c i cos(27r<£>) + a C 2 sm(2Tap) + a c3 cos(47r<£>) + a C 4 sin(47r<£>) =m^ + a c - CS(<p), 

(1) 

where the dot "■" denotes a scalar product. The fourvector CS(y?) is given by the relation: 

CS(y?) = [cos(27r<yj), sin(27r<yj), cos(47r<£>), sin(47r<£>)]. (2) 

The shape of the light curve in the the colour c is then described simply by a fourvector a c , 
a c = [ a ci, a c2, a C 3, a c4 ], consequently by a point in the 4D-space of parameters. The shapes 
of the set of uvbyH p light curves of each of the 82 mCP stars are fully described by the set 
of 5 fourvectors which can be considered 4 matrix Y with 20 elements known with 

different degree of uncertainty. It is rather surprising that there are no two stars with same 
or even resembling set of these vectors. Thus the set of parameters of a given star seem to 
be unique and characteristic. 

3.1 PCA decomposition of light curves of individual stars 

The thorough inspection of uvbyH p light curves of individual stars suggests that these LCs 
are in some extent similar, what entitles us to apply the PCA to them with the aim of 
reducing the number of parameters needed to describe the LCs of a particular star. 

The technique of weighted PCA enables one to search for mutual relations of five fourvec- 
tors of light variations a c in the 4D-space for individual stars. We are looking for the basic 
direction towards which the sum of squared projections is maximum. 

The found orthonormal base vectors d 1; d 2 , d 3 , d 4 determine the course of normalized 
basic LCs fi(<p), ^(v 2 )^ .My?), fM as follows: 



^)=d,-CS(^). 



(3) 
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Figure 2: (a) Two basic phase functions of HD 125248, the first one - fi((p) (full line) and the 
secondary one - /^(v 9 ) (dotted line), (b) The dependence of semiamplitudes A\ c , A 2c and the 
effective amplitude of main and secondary components on wavelength in nm 



We found that a colour light curve m c (ip) can be expressed as a linear combination of at 
most the first three basic functions fj{(f)'- 

k I k \ 

m c {ip) =rn; + J2 A jcfj{ l P) =m;+ ^A, c d, ■ CS(^), k<3. (4) 
i=i \i=i / 

The first basic function, fi(ip) represents common features of the course of all 5 light curves, 
while the second, f2{ ( p), expresses differences from each other. The need of the third function, 
f3{<p), occurs seldom. The accuracy of 0.001 mag of the fitting of the LCs could be achieved 
with only one function f\ in 34 cases, we had to use two functions f\ and f 2 (see Fig. [2h) 
in 47 cases, and we needed three functions fi, f 2 and f% only in the case of HD 62140 (see 
Fig. [3b). 

The semiamplitudes Aj c describe LCs in particular colours quantitatively and give an 
information on mechanisms of stellar variability (see Fig. [2b). They can be determined simply 
as projections of the vector a c into the vector dj of the vector base by the relation: Aj c = a c 
or better by using LSM or a robust regression. We have applied the latter as it suppresses the 
influence of outliers and furnishes us also with the robust values of rrT c for particular colours 
and observers and by the estimate of the uncertainties of the found parameters. 



4 Effective amplitude of a light curve 

The effective amplitude is a robust quantity characterizing the amplitude of a light curve. 



Mikulasek et al.(2007) defined the effective amplitude A e g c of a LC of an mCP star in the 



colour c as 



A eSc = \l8l (m c ((p) - WQ) 2 dip . (5) 
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Figure 3: All the LCs of HD 126515 can be satisfactorily expressed by the linear combination 
of two basic light curves. HD 62140 is the only star where we need three basic light curves. 

The factor 8 is selected so that the effective amplitude of the sine light curve corresponded 
to the amplitude of the observed one. Using relations (T5])and (j3J) we arrive at relations: 
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Here we extended the definition Eq. (jSj) by introducing the negative effective amplitude 
for a LC going in antiphase with the majority of the remaining LCs. The meaning of the 
phrase "to be in antiphase" is specified by the last relation in the equation (jBJ). 

Fig. H] displays that there is no apparent correlation between the mean value of effective 
amplitudes and B-V indices. We argue that mCP stars of various types form more or less 
homogeneous group from the view of their photometric variability. 



4.1 Stars with the largest effective amplitudes 

The median of the set of the effective amplitudes is 0.027 (!) mag. Still some of the mCPs 
have much larger amplitude as documented in the following: the Olsen's star, HD110956B 
has an effective amplitude 0.24 mag in v. The Babckock's star HD 215441 has 0.19 and 0.14 
mag in v and b respectively. HD 51418 has 0.14 and 0.17 mag in H p and y respectively. 
HD 28843, the most variable He-weak star has an effective amplitude 0.15 mag in u. The 
common characteristics of these "record holders" is a single-waved LC (see Figj5]) suggesting 
a simple geometry of their photometric spots. 

4.2 PCA of the dependence of effective amplitudes on wavelength 

The advantage of the PCA is that it offers the possibility to investigate all the sample to be 
studied simultaneously. The light variability of each star can be characterized by a fivevector 
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Figure 4: Dependence of the mean value of the effective amplitude on the average Hipparcos 
B- V. Position of several famous mCP stars are denoted by their HD numbers. We do not see 
any unambiguous relation between the effective amplitudes and B- V indices. 
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Figure 5: The comparison of LCs of two mCP stars with the largest amplitudes. While in the 
case of HD 215441 the amplitudes decrease with the increasing wavelength, right, HD 110956B 
has the largest amplitude in v, and the antiphase occurs in y. These two stars show two basic 
models of the dependence of the amplitude on the wavelength, and represent two different 
mechanisms of variability of the mCP stars. 
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Figure 6: Three most important principle components of the dependence of the effective 
amplitude on the wavelength. The first component occurs in all the sample stars and we 
ascribe it to a bright photometric spot or spots. The second component occurs rather in the 
cooler mCP stars and corresponds to the incidence of a spot extremely dark in the v colour. 
We argue it is the consequence of the diminishing of the Balmer jump in heavily contaminated 
atmospheres. The third component occurs only in a few stars having anomalous appearance 
of the LC in u. 



of the effective amplitudes in all five colours as A cS = [A eSu , A ef[v , A eSb , A e{iHp , A eSy ]. If a LC 
is in antiphase, what occurs approximately in 3.5% of the LCs, the corresponding effective 
amplitude is assigned with the negative sign. 

Then these fivevectors A c g are, using the PCA, expressed by the linear combinations of 
at least as possible number of the main, mutually orthogonal vectors. It turned out that the 
fivevectors A e ^ of all the stars in the sample can be sufficiently accurately expressed as the 
linear combination of three main components (3D-space). Their occurrence and wavelength 
behaviour may reflect some astrophysical significance. Their course with the wavelength is 
displayed on Fig. [6j 

Having the main components defined, we project the vector of effective amplitude of 
an mCP into the 3D basis, thus getting three components b±, 62, &3- They show a specific 
distribution with the B — V index. 



5 Some empirical finding and discussion 

We confirm that the observed light curves of mCP stars can be satisfactorily well represented 
by harmonic polynomials of the second order. Investigating the wavelength dependence of 
the effective amplitude by means of PCA we conclude: 

1. Obviously every of the sample stars has one or more bright spots on its surface whose 
contrast decreases with increasing wavelength. The spots are usually more pronounced 
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Figure 7: The dependence of the projection b\ of the vector of the effective amplitude on 
the first principal component on the B-V index. All the values are positive suggesting the 
photometric spot measured by this quantity looks bright in all colours studied. 
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Figure 8: The dependence of 62 on B — V index. For the hot mCPs the value is close to zero, 
for the cooler ones 62 acquires considerable non-zero values. The corresponding mechanism of 
energy redistribution is most prominent in v colour (see Figj6]), where we see spots in dark. If 
that spot interferes with bright spots culminating near the zero phase it reduces the effective 
amplitude (6 2 > cf. HD 125603). If the dark spot is placed on the opposite side of the star, 
it enlarges the amplitude (6 2 < cf. HD 110956B). 
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Figure 9: The comparison of LCs of two mCP stars with positive hi- In the case of HD 51418 
the amplitude of LC in v is due to dark spot suppressed practically to zero, while in the case 
HD 125630 we even see the v LC to be in antiphase with other LCs. 
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Figure 10: The dependence of 63 on B — V index. For the hot mCPs the value is close 
to zero. For the cooler ones 62 acquires considerable non-zero values and the mechanism 
corresponding to this quantity plays a more significant role. The third principal component 
(see FigJHD acquires large negative value in u implying the spot in this colour is dark. 
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on hotter stars. Nevertheless, occurrence of these spots is not bound to any of the 
peculiarity types. 

2. On cooler mCPs, not necessarily of SrCrRE type, one or more spots looking dark in 
the v colour occur relatively frequently. These spots are indicated by the increasing 
amplitude of the Ci-index in the cool mCPs. We suggest that these spots might be 
related to the decrease of Balmer jump in the spectral energy distribution in regions 
heavily contaminated by overabundant elements. 

3. The spots contrast in the u may demonstrate themselves either as bright or dark, 
whereas as a rule they are bright on hot and dark on cool stars. 

The brightness contrast of the spots probably is due to the uneven distribution of chemical 
elements, horizontal and possibly vertical, and perhaps also to strong global magnetic field. 
These dissimilar atmospheric patches interact differently with the isotropic radiative flux 
coming from the interior. Especially, the "colour" photometric spots may origin due to the 
different redistribution of the flux from the UV to visual region caused by the bound-free 
transitions and line-blocking. Silicon may play a decisive role in the light variability, as it is 
usually overabundant in all subtypes of the mCP stars. 

Consequently, photometric spots are most likely related to the surface distribution of 
chemical elements and, possibly, with the configuration and the intensity of magnetic fields. 
The chemical elements use to be distributed in several more or less irregularly distributed 



spots (cf. Khokhlova et al.(1997)), as well as the magnetic field configuration may differ 
from an elementary dipole structure (IKochukhov 20051) . Thus, a simple pattern of the pho- 
tometric spots cannot be expected, as it was documented e. g. in the case of HD 37776 



(Krticka et al.(2007)) 



The amazing diversity of the forms of the LCs is the result of the complexity of the photo- 
metric spots structure. The study of the relations between the photometric and spectroscopic 
spots on the mCP stars seems to be very promising. 

From this point of view, the mCP stars represent a relatively homogeneous group. How- 
ever, our conclusions are very preliminary mainly for they are based only on the analysis 
of the effective amplitude. We plan also to take into account the shape of a LC harbouring 
information on location of the spots along the surface. 
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